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INTRODUCTION 


Warhead  designs  use  explosives  to  accelerate  components  to  high  velocities.  To 
measure  the  effectiveness  of  an  explosive,  the  cylinder  expansion  test  is  used.  The 
test  measures  the  wall  acceleration  of  a  copper  cylinder  after  the  explosive  inside  has 
detonated.  The  Gurney  constant  is  calculated  from  the  results  and  is  used  to  predict 
the  velocity  in  other  geometry’s.  The  results  of  the  test  are  also  used  in  the  determina¬ 
tion  of  equations  of  state  required  in  computer  modeling  of  warhead  systems. 

In  the  research  and  development  of  new  explosive  molecules  and  formulations, 
quantities  of  materials  are  frequently  limited.  Small  laboratory  batch  processes  would 
have  to  be  repeated  a  number  of  times  in  order  to  produce  the  quantity  required  for  full 
scale  performance  testing,  making  the  test  extremely  expensive.  In  order  to  decrease 
the  cost  and  make  performance  testing,  a  down  selection  criteria  available  earlier  in 
the  development  process,  a  reduced  scale  cylinder  expansion  test  has  been  devel¬ 
oped. 

The  new  test  uses  smaller  cylinders,  3/4  in.  inner  diameter  (I.D.),  1  in.  outter 
diameter  (O.D.),  and  1 0  in.  long,  than  the  conventional  larger  test,  which  is  1  in.  I.D., 

1 .2  in.  O.D.,  and  12  in.  long.  The  volume  within  the  cylinder  and  the  amount  of 
explosive  required  is  reduced  by  53%.  To  directly  compare  the  results  to  the  larger 
test  a  new  Gurney  equation  was  developed.  The  analysis  for  the  first  time  introduces 
variable  velocities  due  to  the  thinning  of  the  wail  during  expansion.  This  allows 
cylinder  tests  of  different  geometry’s  to  be  compared  to  each  other. 


BACKGROUND 

Gurney’s  Method 

R.  W.  Gurney  developed  a  method  for  determining  the  velocity  of  fragments 
accelerated  by  a  detonation.  The  Gurney  method  is  simplistic,  but  the  results  obtained 
correlate  well  with  experiments.  The  model  assumes  a  linear  velocity  distribution 
throughout  the  reacted  gases,  and  a  constant  velocity  in  the  wall.  This  assumption 
creates  a  partitioning  of  the  available  kinetic  energy  between  the  gases  and  the  wail. 
For  an  explosive  driving  a  cylinder 


V=®’ 


(refs.  1  and  2) 


(1) 
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where 


V  =  fragment  velocity 
y2E  =  Gurney  constant  (units  of  velocity) 

M 

Q  =  wall  mass  to  explosive  charge  ratio 

A  useful  engineering  approximation  assumes  the  Gurney  value  to  be  a  constant 
for  the  material.  For  this  reason,  it  is  referred  to  as  the  Gurney  constant.  The  actual 
value  varies  and  is  related  to  many  other  variables,  such  as  the  level  of  expansion, 
brake-up  of  the  wall  releasing  detonation  products,  the  axial  flow  of  the  reaction 
products,  and  the  shock  waves  from  the  detonation. 

The  equation  assumes  that  the  detonation  wave  impacts  the  wall  simultaneously 
along  its  surface.  For  a  cylindrical  geometry  the  detonation  wave  normally  travels 
along  the  axis  of  the  charge.  The  detonation  adds  an  axial  velocity  component.  The 
relationship  between  the  velocity  and  its  normal  component  is  given  by 

V  =  2Dsin(0/2)  (refl)  (2) 


where 


D  =  Detonation  velocity 

9  =  Angle  of  the  tube  wall  to  the  initial  position 

For  many  engineering  problems  all  that  is  required  to  obtain  a  good  estimate  of 
fragment  velocities,  is  to  obtain  a  reasonable  value  for  the  Gurney  constant  from 
testing  or  other  methods.  Equation  2  is  for  a  cylindrical  geometry,  but  with  similar 
assumptions  the  equations  for  many  other  geometries  can  be  obtained.  The  Gurney 
constant  calculated  from  one  geometry  can  be  applied,  with  varying  levels  of  success, 
to  other  geometries. 


PICATiNNY  ARSENAL  CYLINDER  TEST 


In  measuring  the  Gurney  constant,  the  cylindrical  geometry  is  most  often  used.  It 
is  relatively  accurate,  easy  to  manufacture,  and  easy  to  obtain  measurements  from.  In 
the  conventional  test,  as  described  in  reference  3,  the  cylinder  is  made  from  soft 
oxygen  free  copper.  The  wails  do  not  instantaneously  obtain  the  maximum  velocity  as 
assumed  in  the  Gurney  analysis.  These  tubes  allow  for  large  expansions  of  the 
plastically  deformable  copper  before  rupturing.  In  the  conventional  test,  the  tubes  are 
1  in.  I.D.,  1.2  in.  O.D.,  and  12  in.  long.  Other,  usually  larger,  cylinders  have  been  used. 
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The  Picatinny  Arsenal  cylinder  expansion  test  uses  a  slightly  different  geometry 
(fig.  1).  The  slightly  smaller  tubes  are  0.75  in.  I.D.,  approximately  1  in.  O.D.,  and  10  in. 
long.  The  smaller  geometry  uses  a  little  less  than  half  of  the  explosive  required  for  the 
larger  test.  This  is  a  significant  savings  for  research  materials  that  are  made  in  very 
small  quantities.  The  walls  are  thicker  than  the  conventional  test.  This  allows  for 
further  expansions  before  breakup  and  changes  the  charge  to  mass  ratio  (C/M)  of  the 
test.  A  typical  Picatinny  Arsenal  test  will  have  a  C/M  of  about  0.25,  which  would  be 
0.442  for  the  conventional  cylinder  expansion  test  for  an  explosive  of  the  same 
density. 


Figure  1 

Picatinny  Arsenal  cylinder 

The  wall  and  detonation  velocities  are  obtained  in  the  same  manner  as  for  the 
conventional  tube.  A  streak  camera  is  used  to  record  the  expansion  of  the  cylinder 
walls.  The  camera  is  typically  run  at  a  writing  rate  of  1  mm/ps,  allowing  the  film  image 
to  fit  on  the  screen  of  the  motion  analyzer.  The  image  of  the  copper  cylinder  is  located 
in  the  center  of  the  field  of  view,  allowing  for  the  analysis  to  be  conducted  on  both 
sides  of  the  tube.  The  slit  of  the  camera  is  located  approximately  6.5  in.  from  the  top  of 
the  tube.  This  places  the  slit  far  enough  from  the  ends  to  avoid  rarefactions  from 
effecting  the  results.  If  the  explosive  fill  consists  of  individual  pellets,  the  slit  is  moved 
until  it  is  located  away  from  the  edges  of  any  of  the  pellets.  This  avoids  any  detonation 
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transfer  problenns  and  any  air  in  this  region  from  adversely  effecting  the  results.  To 
assure  the  detonation  is  well  established  before  entering  the  cylinder,  2  to  3  in.  of  run 
up  are  used.  This  long  run  up  replaces  the  expensive  plain  wave  generator  that  is 
sometimes  used. 

The  tube  is  back  lighted  with  an  argon  bomb  (fig.  2).  For  this,  an  argon  filled  toy 
balloon  is  used  (ref  4)  with  an  explosive  charge  taped  to  the  back.  The  balloon 
provides  about  50  ps  of  usable  light  and  is  cheap  and  easy  to  use.  This  also  reduces 
the  cost  as  compared  to  the  more  expensive  argon  bombs  typically  used. 


Figure  2 

Experimental  setup 

The  detonation  velocity  is  obtained  using  closure  pins,  placed  at  1  and  2  in.  from 
either  end  of  the  tube  (fig.  3  and  discussions  with  I.  Akst).  In  this  system,  the  copper 
tube  is  grounded.  Ordinance  tape  is  wrapped  once  around  the  tube  and  copper  wire 
of  0.004  in.  diameter  with  cleaned  ends  is  placed  on  the  tape  extending  3/16  in.  over 
the  side.  The  wire  is  then  taped  in  place.  This  puts  the  copper  wire  laying  3/16  in. 
from,  but  not  touching,  the  copper  tube.  The  pins  are  then  charged.  Contact  between 
the  pins  and  the  copper  tube  during  the  test  is  measured  with  a  capacative  discharge 
system.  An  oscilloscope,  recording  at  5  ris  per  point,  measures  the  shorting  of  the 
pins. 
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Figure  3 

Closure  pin  locations 


Analysis 


The  streak  photograph  records  the  position  of  the  wall  with  time  (fig.  4).  By 
scaling  the  photograph,  the  position  time  history  in  full  scale  is  obtained.  A  Vangaurd 
brand  motion  analyzer  is  used  to  obtain  the  data,  which  is  transferred  into  an  IBM 
compatible  personnel  computer.  Two  programs  are  used,  the  first  scales  and  saves 
the  data  in  a  file  and  the  second  calculates  the  Gurney  values.  This  allows  for  any 
future  modifications  in  the  analysis  without  re-entering  the  data,  a  time  consuming 
process.  Generally,  about  100  to  200  readings  are  obtained  on  the  cylinders  expan- 
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Figure  4 

Typical  test  results 

The  second  program  analyzes  the  scaled  data.  This  program  performs  a  least- 
mean-square  that  fits  the  data  to  a  fourth  order  polynomial.  From  the  differential  of  the 
fitted  equation,  the  velocity  can  be  obtained  at  any  point  in  time  up  to  break-up. 

The  velocity  varies  in  the  cylinders  during  the  expansion.  The  wall  thins,  giving 
the  inner  wall  a  slightly  higher  velocity  than  the  outer  wall.  The  Picatinny  Arsenal 
cylinder  test,  with  its  heavier  wall,  has  a  larger  variation  of  velocities  than  the  conven¬ 
tional  tube.  A  modified  Gurney  equation  was  developed  that  accounts  for  these 
variations,  and  is  discussed  later  (Gurney  Relationship  for  Thick  Cylinders). 

The  program  determines  the  velocity,  and  therefore  the  Gurney  constant,  at  two 
and  seven  volume  expansions  of  the  inside  diameter.  These  are  the  same  expan¬ 
sions  as  used  in  conventional  expansion  tests.  The  Picatinny  Arsenal  cylinder  has  not 
broken  up  as  yet  at  much  larger  expansions,  obtaining  at  least  11.7  volume  expan¬ 
sions  of  the  inside  diameter;  but  the  Gurney  values  obtained  are  only  slightly  higher 
than  those  obtained  at  the  seven  area  expansions.  The  seven  volume  expansions  of 
the  outside  diameter  are  typically  used  since  this  is  the  volume  expansion  reliably 
obtained  before  the  conventional  test  breaks  up.  The  break  up  is  readily  evident  as 
the  gases  jetting  from  the  cracks  make  the  cylinder  appear  to  be  rapidly  accelerating 
after  a  long  period  of  steady  velocity.  The  typical  output  from  the  program  is  in  the 
appendix. 
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VELOCITY  RELATIONSHIPS  OF  CYLINDER  EXPANSION  TESTS 
Purpose  for  Relationships 

For  the  cylinder  expansion  tests,  the  velocity  measurements  are  obtained  on  the 
outside  diameter  of  the  tube.  Because  of  the  thinning  of  the  tube  during  expansion, 
the  inside  wall  has  a  higher  velocity  than  the  outside.  Also,  the  inside  area  expands  at 
a  rate  much  higher  than  the  area  measured  on  the  outside  wall.  For  this  reason,  some 
examination  of  the  area,  velocity,  and  other  relationships  during  the  expansion  are 
useful. 

Momentum  Average  Velocity 

For  a  cylinder,  the  mass  M  remains  constant  during  the  expansion 


M  =  pr^7i-prf7t  =  p7t(r^-tf) 


Change  in  mass  with  time  is  zero 


riM  dr. 

^  =  0  =  p7i  2  r^^  -  r. 
dt  ^  L  o  dt  '  dt 


Then 


dr  dr, 

r  — -  -  r. — - 

°dt  ‘dt 


V  r.  V 

— 2-  =  -L  and  for  any  radius  within  the  cylinder  -^  =  -L 
V.  r„  .  V  r 


The  average  velocity  of  the  copper  wall  must  be  calculated  based  upon  the  mass 
or  the  energy  of  the  wall.  The  first  computation  here  uses  the  momentum  average 
velocity  as  compared  to  the  velocity  measured  on  the  outside  wall. 

For  the  momentum  I 


Vp27r  rdr 
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Using  the  equation  7  relating  the  velocities  and  radii 


pv 


27i:rdr 


The  momentum  average  velocity  is  then 

V  =  ^0^0  (*^0  '  ^i)  [Momentum 

(r^-rf)2ip  I  Weight  . 

■  (^-f) 


(8) 


(9) 


From  the  constraint  of  constant  mass 

1^  -  r^  =  constant  for  the  system 

(10) 

=  (^-f)|  INITIAL 

Then 

^  =  I  initial)^  (11) 

Then  Vm  becomes 

V  -  In I  initial)^) 

I  initial  (12) 

The  ratio  of  the  momentum  average  velocity  to  the  outside  wall  velocity  is  shown 
in  the  following  graph.  The  initial  conditions  are  for  the  Picatinny  Arsenal  cylinder  test. 
As  the  cylinder  expands  and  the  walls  become  thinner,  the  momentum  average 
velocity  and  the  outside  wail  velocity  match  closely. 
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Figure  5 

Ratio  of  momentum  average  velocity  to  outside  wall  velocity 

This  yields  the  momentum  average  velocity  from  the  velocity  of  the  outside  wall 
and  the  outside  radius.  Note  that  for  expansions  in  the  areas  of  interest,  the  momen¬ 
tum  average  velocity  and  the  outside  wall  velocity  are  essentially  the  same. 

Kinetic  Energy  Average  Velocity 

Looking  at  kinetic  energy 


KE  =  1/2mV^ 


V  pTi  rdr 


forV  = 


(13) 


ro .  ,2 


Jr-, 


vV 

o  o 

r^ 


pTT  rdr 


Solving  and  using  a  equation  10 


KE  =  pvfr^  In 


(14) 


Which  yields  the  average  kinetic  energy  from  the  velocity  measured  on  the  outside  of 
the  radius  and  the  outside  diameter. 
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The  kinetic  energy  average  velocity 


V„ 


2  KE 


Mass 


2  V^r^  In  I 
0  0 

1  ^0  1 

|2 

initial))) 

INITIAL  ) 

(15) 


The  kinetic  average  velocity  divided  by  the  outside  wall  velocity  is  then 


1 

>  1 

2  1^  In  ^ 
0 

\ 

r  \ 

0 

('^-(('^■'f)!  initial))^) 

<1 
O  1 

1 

((■^  -  -f)  INITIAL  ) 

(16) 

The  following  diagram  shows  the  relationship  between  the  kinetic  average 
velocity  and  the  measured  outside  wall  velocity.  The  same  trends  are  shown  as  for  the 
momentum  average  velocity,  the  ratio  nears  one  as  the  walls  thin  due  to  expansion. 


0  2  4  6 

radius/  initial  radius 

Figure  6 

Ratio  of  kinetic  energy  average  velocity  to  outside  wall  velocity 


Measurement  Errors 


A  rigorous  analysis  of  the  errors  measurements  have  not  been  conducted.  The 
analysis  here  determines  the  general  relationships  between  the  variables. 
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For  the  equation  relating  velocity  and  charge  to  mass  ratio,  C/M 


1 


It  can  be  seen  that  for  smaller  M/C  ratios,  the  errors  in  determining  the  errors  are 
not  as  significant  if  the  M/C  ratio  is  small.  The  errors  in  measuring  the  velocity  have  a 
much  greater  effect,  linearly  changing  the  calculated  value  of  the  Gurney  constant. 

Error  From  Not  Accounting  for  the  Angle  of  Expansion 

From  the  relationship  between  angle,  wall,  and  detonation  velocities 
V  =  2Dsin(e/2) 


Where  V  is  obtained  from  the  Gurneys  relationship 

The  measurement  used  in  the  test  is  the  radial  position,  so  only  radial  velocities  are 
obtained. 


V  ,  =  V  cos  F 

measured  \2 


Using  trigonometric  identities 


V„sasu,ed  =  °  S'"(e| 


Then 


=  2D  sin  ||-|  cos  ||- 


6  =  sin'^  [^measured 

\  D 


and 


V  = 


V 


measured 


.  -1 
COS  isin 


V 


measured  I 


D  / 


(17) 


(18) 


(19) 


(20) 


For  most  problems,  this  accounts  for  a  small  percent  of  error  between  the  actual  and 
the  measured  velocities.  This  correction  factor  can  then  be  ignored. 
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Gurney  Relationship  for  Thick  Cylinders 

Previous  work  on  the  Gurney  constant  assumes  that  the  wall  thickness  of  the 
cylinder  remains  constant.  This  is  not  correct  as  the  wall  must  become  thinner  as  the 
tube  expands.  To  correct  for  the  slightly  larger  inside  wall  velocity,  some  researchers 
have  applied  a  momentum  correction  to  the  values  obtained.  It  is  possible  to  obtain  an 
expression  that  does  not  require  such  a  fix,  by  applying  the  conservation  of  mass. 

Velocity  Distribution  of  Gasses  and  Wall 


By  the  conservation  of  mass 
V.r 

V  =  —  V  =  Gas  velocity  at  radius  r 

gas  f.  gas  ’ 

V  r 

^waii  “  ^waii  “  velocity  at  radius  r 

Vo  =  Outside  wall  velocity 
Vj  =  Inside  wall  velocity 
ro  =  Outside  wall  radius 
n  =  Inside  wall  radius 
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By  the  conservation  of  energy 


CE  =  1/2mV  .,^  +  1/2mV„^/ 

wall  gas 


v  =  f^ 

L  \  O  I  !  J 


(21) 
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Weight  of  the  wall,  ^  =  Pm  ^ 
Weight  of  the  explosive,  C  =  ti 


The  M/C  ratio  can  be  obtained  from  the  equation  to  have  it  resemble  the  typical 
form  in  equation  1,  and  be  based  upon  the  outside  radius  position,  the  measured 
quantity  in  testing. 

From  the  conservation  of  mass 

^  =  initial)^ 


Then 


V= 


(tl'il^klNITIAL 


pM 

0^0  1  ^ 

INITIAL  ^ 

Inj 

/  ^ 
c 

\ 

) 

1  (^^)l 

INITIAL  1 

1 

)2 

INITIAL  /  / 

+  1/2 


(22) 


This  equation  can  be  used  to  determine  the  Gurney  constant  for  cylinders  with 
thick  walls.  It  also  can  be  used  to  compare  the  test  results  with  different  cylindrical 
geometry's.  The  equation  was  manipulated  to  be  dependent  upon  the  outside  radius 
and  the  initial  conditions  of  the  tube.  This  allows  for  easy  analysis  with  the  measured 
variables.  A  comparison  for  the  Picatinny  Arsenal  cylinder  test  in  given  in  figure  7. 
The  effects  of  this  equation  appear  more  in  the  early  expansion  of  the  tube. 
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Ratio  of  Modified 
Gurney  Velocity 
to  Conventional 


Ratio  of  inside  cross  Sectional 
Area  to  Initial  Area 


Figure  7 

Ratio  of  modified  gurney  to  gurney  constant 


SUMMARY 

A  smaller,  relatively  heavier  walled  copper  cylinder  has  been  incorporated  into 
the  cylinder  expansion  test.  Modifications  to  the  Gurney  Equation  that  account  for  the 
stretching  and  thinning  of  the  copper  cylinder  have  been  developed.  The  new 
equation  allows  the  use  of  the  smaller  cylinders,  3/4  in.  inner  diameter  (I.D.).,  1  in. 
outter  diameter  (O.D.).,  and  10  in,  long  to  the  conventional  larger  test,  1  in.  I.D.,  1.2  in. 
O.D.,  and  12  in.  long.  This  reduces  the  quantity  of  explosives  required  for  performance 
testing  by  more  than  half. 
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APPENDIX 

INPUT  AND  ANALYSIS  PROGRAMS 
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Input  Program 

This  program  retrieves  the  raw  data  from  the  serial  port,  collects 
scales  and  other  test  data,  and  saves  the  data  in  scaled  format. 

#include  <stdio.h> 

#include  <math.h> 
double  x[300],y[300]; 
char  header[8] , infile  [  1 2]  ,outfile[  1 2] ; 
main()  { 

FILE  *data_in; 

FILE  *data_out; 
int  rawcount,i; 

char  xstring[7],ystring[7],tmp[7]; 

double  filmspeed,od,id,moverc, timescale, radiuscale,rho,DefVel; 
rawcount=0; 

printfC'Key  in  the  input  file  name  "); 
scanf("%s",infile); 

printf("\nKey  in  the  output  filename  "); 
scanf("%s",outfile); 
data_in=fopen(infile,  "r  ") ; 
for(i=0;  i<299 ;i++) 

{if(feof(data_in)!=  0)  break; 
fscanf(data_in,"%s%s%s",header,xstring,ystring); 
x[i]=atof(xstring);  y[i]=atof(y string); 

}  fclose(data_in); 
rawcount=i-l; 

printf("Number  of  Data  Points  Read  %d\n",rawcount); 
printf('Nn\n\n\nBldg.  #1612  Run  Number  3-%s",&(header[5])); 
printf('^n\n\n\nEnter  The  Following  Variables\n\n"); 
printfC'The  M/C  Ratio  ");scanf("%s",tmp); 
moverc=atof(tmp); 

printfCNnThe  O.D.  of  the  Tube(inches)  ");scanf("%s",tmp); 
od=atof(tmp); 

printf('NnThe  LD.  of  the  Tube(inches)  ");scanf("%s",tmp); 
id=atof(tmp); 

printf("\nThe  Film  Mirror  Speed(R.P.S.)  ");scanf("%s",tmp); 


filmspeed=atof(tmp); 

printf('\nThe  Explosive  Density  (gm/cc)  ");scanf("%s",tmp); 
rho=atof(tmp); 

printfC  ViThe  Detonation  Velocity  (mm/us)  ");scanf("%s",tmp); 
DetVel=atof(tmp); 

timescale=fabs(((4*4.75)/(x[0]-x[l]))*(250/filmspeed)); 

radiuscale=fabs(od*25.4/(y[2]-y[3])); 

for(i=4;  i<ra  wcount;  i++) 

{x[i-4]=x[i]*timescale;  y[i-4]=y[i]*radiuscale; } 
rawcount=rawcount-4; 
data_out=fopen(outfile,"w"); 

fprintf(data_out,"3-%s  Run  Nuniber\n\n'',&header[5]); 

fprintf(data_out,"%f  Mass/Charge  Ratio\n",moverc); 

fprintf(data_out,"%f  Outside  Diameter  nim\n",od*25. 4); 

fprintf(data_out,"%f  Inside  Diameter  mm\n",id*25.4); 

fprintf(data_out,"%f  Explosive  Density  gm/cc\n",rho); 

fprintf(data_out,"%f  Detonation  Velocity  mm/us\n",DetVel); 

fprintf(data_out,  "%d  Total  Number  of  Data 

Point  s\n '  ',ra  wcount) ; 

fprintf(data_out,'\n\nTime(us)  Radius(mm)\n\n"); 
for(i=0;i<rawcount;i-f+)  fprintf(data_out,"%f  %fm",x[i],y[i]); 
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Analysis  Program 

This  program  analysis’s  the  data  retrived  by  the  previous  program. 

#include  <stdio.h> 

#include  <math.h> 

#include  <stdlib.h> 

#defme  PI  3.14159 

double  x[300],y[300],arry[5][6]; 

double  position[300],vel[300],accel[300],gur[300],volf[300]; 

char  runnum[6],infile[  1 2],outfile[  1 2]  ,outreport[  1 2] ; 

double  bisect_time  (double, double[5]  [6], double); 

double  velocity_of_wall(double,double[5][6]); 

double  position_of_wall(double,double[5][6]); 

double  modified_gumey(double, double, double, double, double); 

main()  { 

FILE  *data_in; 

FILE  *data_out; 

FILE  *report_out; 

int  count,!, n,m,nc,mc,npm,j,k; 

char  *trash,*tmp; 

double 

od,id,moverc,rho,DetVel,err,deter,t,volexp,Ydesired,bitime,v,mgur; 

printfC'Key  in  the  input  file  name  ");  scanf("%s",infile); 
data_in=fopen(infile, '  'r' ') ; 

p^intf('^nKey  in  the  output  file  name  ");  scanf("%s",outfile); 
printf("\nKey  in  the  report  file  name  ");  scanf("%s",outreport); 
report_out=fopen(outreport,  "w  ") ; 
fscanf(data_in,"%s%s%s",runnum,trash,trash); 
printf('\n  Bldg.  1612  Run  Number  %s",runnum); 
fprintf(report_out,"Bldg.  1612  Run  Number  %s",runnum); 
fprintf(report_out,'NnFrom  data  file  %s",outfile); 
fscanf(data_in,"%s%s%s",tmp,trash,trash); 
moverc=atof(tmp); 

fscanf(data_in,"%s%s%s%s",tmp,trash,trash,trash); 

od=atof(tmp); 

fscanf(data_in,"%s%s%s%s",tmp,trash,trash,trash); 

id==atof(tmp); 

fscanf(data_in,"%s%s%s%s",tmp,trash,trash,trash); 
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rho=atof(tmp); 

fscanf(data_in,"%s%s%s%s",tmp, trash, trash, trash); 
DetVel=atof(tmp); 

fscanf(data_in,"%s%s%s%s%s%s",tmp,trash,trash,trash,trash, trash); 
count=atoi(tmp); 

fscaiif(data_in,"%s%s", trash, trash); 
for(i=0;i<couiit;i-f +)  { 

fscanf(data_m, '  '%s  ",tmp) ; 

x[i]=atof(tmp); 
fscanf(data_in,  "%s '  ',tmp) ; 
y[i]=atof(tmp); } 
fclose(data_m); 

/*initialize  array*/ 
for  (n=0;n<5;n++) 

{for(m=0;m<6;m++) 

arry[n][m]=0;} 

/♦Heart  of  the  Array*/ 
for  (n=0;n<5;n++) 

(for  (m=0;m<5;nn-+) 

{for(i=0;i<count;i++) 

if(m+n==0)  arry[n][m]=arry[n][m]+l; 

else  {arry [n]  [m]=arry [n]  [m]+pow(x[i] ,n+m); )  }  } 
/*Vector  of  the  Array*/ 
for  (i=0;i<count;i++) 

{for(n=0;n<5  ;n++) 
if(n==0)  arry[n][m]=arry[n][5]+y[i]; 

else  {arry[n][5]=arry[n][5]+y[i]*pow(x[i],n);}} 
printf("\n"); 

/*  simeq  */ 

nc=4; 

mc=l; 

npm=nc+mc; 

err=.0000000001; 

deter=l; 

for(i=0;i<=nc;i-f+) 

{deter=deter*arry  [i]  [i] ; 

for(j=i+ 1 ; j<=npm;j++)  arry [i]  [j]=arry [i]  Ij]/arry [i]  [i] ; 

arry[i][i]=l; 

for(j=0;j<=nc;j++) 
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{if(i!=j)  {if(arry[j][i]!=0) 

{for(k=i+ 1  ;k<=npm;k++) 

arrylj]  [k]=aiTy  [j]  [kj-arry  Ij]  [i]*arry  [i]  [k] ; 
any|j][i]=0;}>}> 

fprmtf(report_out,'Nn\n  Array  Solution  to  Least  Mean  Square  Fit\n"); 
for(i=0;i<=nc;i++) 

{fprintf(report_out,  "\n ") ; 
for(j=0;j<=npm;j++) 
fprintf(report_out,"  %f' ,arry[i][j]); } 

/*  Gurney  Calculations  */ 
data_out=fopen(outfile, ' '  w ' ; 
for(i=0;i<count;i++) 

{t=x[i]; 

position  [i]=position_of_wall  (t,arry); 
vel[i]=velocity_of_wall  (t,arry); 
accel[i]=2*arry[2][5]+6*arry[3][5]*t+12*arry[4][5]*t*t; 
gur[i]=vel[i]/(pow(moverc+.5,-.5)); 
volf[i]=(pow(y[i]*2+od,2)-pow(od,2)+pow(id,2))/pow(id,2); 
mgur=modified_gumey(vel[i],od,id,moverc,y[i]); 
fprintf(data_out,"%f  %f  %f  %f  %f  %f  %f  %f\n" 
»y  [i]  ,position[i] ,  vel[i]  ,accel[i]  ,gur  [i]  ,x  [i] ,  volf  [i]  ,mgur) ; } 
fclose(data_out); 
printfCNn"); 
for  (i=2;i<=7;i++) 

{volexp=i*PI*pow(id/2,2)+(PI*(pow(od/2,2)-pow(id/2,2))); 

Ydesired=pow(volexp/PI,,5)-od/2; 

bitime=bisect_time  (Y desired, any, x [count- 1  ]); 

printf('\n"); 

fprintf(report_out, '  Nn") ; 

printf("\nFor  %d  Volumes  Expansion  Time  =  %f  us",i,bitime); 
fprintf(report_out,'\nFor  %d  Volumes  Expansion  Time=  %f 
",i,bitime); 

v=velocity_of_wall(bitime,arry); 
printf('NnWith  a  Wall  Velocity  of  %f ',v); 
fprintf(report_out,"\nWith  a  Wall  Velocity  of  %f  mm/us",v); 
printf('ViAnd  a  Gurney  of  %f  mm/us", v/(pow(moverc+.5,-.5))); 
fprintf(report_out,'NnAnd  a  Gurney  of  %f 
mm/us",v/(pow(moverc-h.5,-.5))); 

mgur=modified_gumey(v,od,id,moverc,Ydesired); 
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printf('\nAnd  a  Modified  Gurney  of  %f  min/us",mgur); 
fprintf(report_out,'\nAnd  a  Modified  Gurney  of  %f  mm/us  ",mgur) ;  > 
fclose(report_out); 

> 

/* _ _ _ _ _ _ _ */ 

double  bisect_time  (double  Ydesired,  double  any  [5]  [6], double  Tupper) 
{double  Ygeuss,Tlower,Tmid; 

Tlower=0; 

while(floor(Ygeuss*  1 00000)  !=floor(  Ydesired*  100000)) 
{Tmid=(Tupper-fTlower)/2; 

Y  geuss=arry  [0]  [5]+arry  [  1  ]  [5]  *Tmid+arry  [2]  [5]  *pow(Tmid,2)-i-arry  [3 
][5]*pow(Tmid,3)+arry[4][5]*pow(Tmid,4); 
if(Ygeuss  >  Ydesired)  Tupper=Tmid; 
if(Ygeuss  <  Ydesired)  Tlower=Tmid; } 
return  Tmid; } 

/* _ _ _ */ 

double  position_of_wall  (double  t,  double  arry[5][6]) 

{double  position; 

position=arry  [0]  [SJ+arry  [  1  ]  [5]  *t+arry  [21  [5]  *t*t-}-arry  [3]  [51  *t*t*t+arry  [4]  [51 

*t*t*t^t* 

return  position; } 

/* _ _ _ _ _ */ 

double  velocity_of_wall  (double  t,  double  any  [5]  [6]) 

{double  velocity; 

velocity =arry[  1  ]  [51+2*arry  [2]  [5]  *t+3  *arry  [3]  [51  *t*t-f  4*arry  [4]  [5]  *t*t*t; 
return  velocity; } 

/* _ _ _ */ 

double  modified  _gumey  (double  v,  double  od,  double  id,  double 
moverc, double  rout) 

{double  gamma; 
double  mgur; 
double  top; 
double  bottom; 
rout=rout+od/2; 

gamma=pow(od/2,2)-pow(id/2,2); 

top=pow(pow(rout*rout-gamma,.5)/rout,2); 

bottom=2*moverc*((rout*rout- 

gamma)/gamma)*log(rout/(pow(rout*rout-gamma,.5)))+.5; 
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mgur=v/po  w((top/bottom),.5); 
return  mgur; } 

/* - - - - .*/ 
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